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Abstract. We consider a one-dimensional diffusion process with coefficients 
that are periodic outside of a finite 'interface region'. The question investigated 
in this article is the limiting long time / large scale behaviour of such a process 
under diffusive rescaling. Our main result is that it converges weakly to a 
rescaled version of skew Brownian motion, with parameters that can be given 
explicitly in terms of the coefficients of the original diffusion. 

Our method of proof relies on the framework provided by Freidlin and 
Wentzell |FW93I for diffusion processes on a graph in order to identify the 
generator of the limiting process. The graph in question consists of one vertex 
representing the interface region and two infinite segments corresponding to 
the regions on either side. 



1. Introduction 

Consider a diffusion process in M.'^ of the type 

dX{t) = b{X)dt + dB{t) , X{0)eR'^, (1.1) 

where B is a d-dimensional Wiener process and 6 : M'' ^ W'' is periodic and smooth, 
and define the diffusively rescaled process X^{t) — eX{t/e^). If h is periodic and 
satisfies a natural centering condition, then it is well-known that X^ converges in 
law as e — to a Wiener process with a 'diffusion tensor' that can be expressed in 
terms of the solution to a suitable Poisson equation, see for example the monographs 
|ijLP78[IP508] . 

Similar types of homogenization results still hold true if b is not exactly periodic, 
but of the form b{X) — b{X^eX), for some smooth function b that is periodic 
in its first argument. In other words, b consists of a slowly varying component, 
modulated by fast oscillations. In this case, the limiting process is not a Brownian 
motion anymore, but can be an arbitrary diffusion, whose coefficients can again be 
obtained by a suitable averaging procedure |BMP05) . The aim of this article is to 
consider a somewhat different situation where there is an abrupt change from one 
type of periodic behaviour to another, separated by an interface of size order one 
in the original 'microscopic' scale. 

To the best of our knowledge, this situation has not been considered before, 
although a similar problem was studied in |ACP03| . In order to keep calculations 
simple, we restrict ourselves here to the one-dimensional situation. Building on this 
analysis, we are able to address the multidimensional case in iHMlOj . Restrict mg 
ourselves to the one-dimensional case considerably simplifies the analysis due to the 
following two facts: 

• Any one-dimensional diffusion is reversible, so that its invariant measure 
can be given explicitly. 
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• The 'interface' is a zero-dimensional object, so that it cannot exhibit any 
internal structure in the limit. 

Before we give a more detailed description of our results, let us try to 'guess' what 
any limiting process should look like, if it exists. Away from the interface, we 
can apply the existing results on periodic homogenization, as in |BLP78[ IPS08) . We 
can therefore compute diffusion coefficients C± such that X° is expected to behave 
like C+W{t) whenever X'^ > and like C-W{t) whenever X'^ < 0, for some Wiener 
process W. One possible way of constructing a Markov process with this property 
is to take = G{W{t)), where G: M ^ M is given by 

C+x if a; > 0, 
C_x otherwise. 



In turns out that processes of this form do not describe all the possible limiting 
processes that one can get in the presence of an interface. The reason why this 
is so can be seen by comparing the invariant measure of G(W{t)) to the invariant 
measure of X"^ . Since the invariant measure for W is Lebesgue measure (or multiples 
thereof), the invariant measure for G{W{t)) is given by 

( X+dx if a; > 0, . A+ C_ 

^'(d^) = i X^dx otherwise, " = " • (1-2) 



A- 



On the other hand, if we denote the invariant measure for X^ by /i^ = p^{x)dx^ 
then will typically look as follows: 



p-{x) 




It follows that one does indeed have ji' 
but the ratio A 



(1.2) 



./A_ 



/i" as e — > 0, where is of the type 
depends on the behaviour of &, not only away from 



the interface, but also at the interface. This can be understood as the process X^ 
picking up an additional drift, proportional to the local time spent at 0, that skews 
the proportion of time spent on either side of the interface. A Markov process with 
these properties can be constructed by applying the function G to a skew-Brownian 
motion (see for example |Lej06| ) with parameter p for a suitable value of p. 

An intuitive way of constructing this process goes as follows. First, draw the 
zeroes of a standard Wiener process on the real line. These form a Cantor set 
that partitions the line into countably many disjoint open intervals. Order them by 
decreasing length and denote by /„ the length of the nth interval. For each n > 0, 
toss an independent biased coin and draw an independent Brownian excursion. If 
the coin comes up heads (with probability p), fill the interval with the Brownian 
excursion, scaled horizontally by /„ and vertically by G+y'^- Otherwise (with 
probability 1 — p) , fill the interval with the Brownian excursion, scaled horizontally 
by /„ and vertically by —C-\/Tn. One can check that the invariant measure for 



this process is given by (1.2 1, but with 

A_ (1 



pG- 



(1.3) 
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We denote the corresponding process by Bc±,p{t). 

This should ahuost be sufficient to guess the main result of this article. To fix 



notations, we consider the process X{t) as in ( 1.1 1 and its rescaled version , and 
we assume that the drift function h is smooth and periodic away from an 'interface' 
region [—77, 77]. More precisely, we assume that there exist smooth periodic functions 
6i : M — M, z G {+, — }, such that hi{x + 1) = bi{x) and such that h{x) = b+{x — rf) 
for X > rj and b{x) — b-{x + rf) for x < —77. Additionally, we assume that the 
functions bi satisfy the centering condition 

»i 

bi{x) dx = . 

We also set V{x) = b{x) dx for x e M, so that exp{2V {x)) dx is invariant for X, 
and similarly for Vi. Denote by Ci the effective diffusion coefficients for the periodic 
homogenization problems corresponding to b^ (see equation (3.7) below or [PS08| 
for a more explicit expression). Define furthermore A± by 

A+ = / exp{2V{x))dx , ^- = cxp{2V(x)) dx , (1.4) 

J 7] J 



and let p e (0, 1) be the unique solution to (1.3). With all these notations at hand, 
we have: 

Theorem 1.1. For any t > 0, the law of X^ converges weakly to the law of Bc±,p 
in the space C{[0,t],R) . 

Remark 1.2. In order to keep notations simple, we have assumed that the diffusion 
coefficient of X is constant and equal to 1. The case of a non-constant, but smooth 
and uniformly elliptic diffusion coefficient can be treated in exactly the same way, 
noting that it reduces to the case treated here after a time change that can easily 
be controlled. 

Remark 1.3. A natural extension of the results presented in this article is case of a 
random potential, along the lines of the situation first considered in |PV82j . While 
the heuristic argument presented in the introduction still applies, the method of 
proof considered here does not seem to apply readily. 

An alternative method would be to consider an injective harmonic function TJ^ 
for X^, so that H^{X'^{t)) is a martingale for every e. In the situation at hand, one 
can take 



H,{x)^ r e^'^^y'^Uy 
Jo 



One might hope that in this case it is possible to show that the process H^{t) — 
H^{X^{t)) then converges as e — > to the martingale H{t) given by 



H{t)= f A{H{s))dB{s) , 
Jo 



where A{H) takes one constant value for > and a different constant value for 
H < (see also Section |4] below) . This type of approach has been successfully 
applied to a number of multiscale problems, including homogenisation on fractals, 
see for example |Koz93l [Ziu95l IKK961 IHK981 IOwh03l IKum04j . 

The proof of the weak convergence of the probability distributions on C[0,oo) 
associated to X| presented in this article will rely heavily on the 1993 paper by 
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Freidlin and Wentzell [FW93| . in which the authors consider a 'fast' Hamihonian 
system perturbed by a 'slow' diffusion. Theorems 2.1 and 4.1 from |FW93j provide 
a general framework for proving first the tightness and then the convergence of a 
family of probability distributions on C[0, oo). 

We will start by showing tightness of our family of processes in Section [2j Once 
tightness is established, we show in Section [s] that every limiting process solves 
the martingale problem associated to a certain generator. Finally, we show in 
Section [4] that this martingale problem has a unique solution which is precisely the 
rescaled skew-Brownian motion, thus concluding the proof. 

Acknowledgements. We would like to thank the anonymous referees for carefully read- 
ing a preliminary version of this article and encouraging us to shorten some of our argu- 
ments considerably. The research of MH was supported by an EPSRC Advanced Research 
Fellowship (grant number EP/D071593/1) and a Royal Society Wolfson Research Merit 
Award. 



2. Proving tightness 

The main result of this section is the following: 

Theorem 2.1. The family of probability measures on C([0, oo), M) given by the laws 
of for e e (0, 1] is tight. 

Proof. Given that the initial condition is kept fixed at one single point across the 
entire family of laws, tightness follows from uniformity in the modulus of continuity 
over £ e (0, 1]. 

For X G E and p > 0, denote by the first exit time of the canonical process 
from the interval \x — p,x + p\. We also denote by f'x,e the law of X^. It then follows 
immediately from the proof of |SV791 Theorem 1.4.6] that a sufficient criterion for 
tightness is that, for every p > 0, there exists a constant Ap such that the bound 

Px,.(t;<^)<m, , (2.1) 

holds uniformly over all x G M, e € (0, 1], and (5 G M. Before we proceed, we note 
the following two crucial facts: 

(1) It follows from the periodic case |BLP781 Section 3.4] that (2.1) holds uni- 
formly for X ^ {—£ri ^ p, ST] + p). We denote the corresponding constants 

by^i. 

(2) A standard martingale argument as in |SV79[ Section 1.4] shows that, for 



every eg > 0, the bound (2.1 ) holds uniformly over all a; G M, provided that 
we restrict ourselves to e G [eq, 1]. For the sequel of the proof it will be 
convenient to make the choice Eq = p/ (Arj) and we denote the corresponding 
constants by A-^. 

Combining these two facts, we see that it remains to find a family of constants 
such that (2.1) holds for x G [—67] — p,er] + p) and for e < p/{Arj). At this stage 



we note that since p is greater than twice the width of the interface, for every a; G M 
there exist two points x± with the following two properties: 

(1) The process started at x has to hit either x+ or a;_ before it can reach the 
boundary of the interval [a; — p, x -f p] . 

(2) The intervals I± — [x± — f,i± -I- |] satisfy I± n (— 776, 775) — and I± C 
[x- p,x + p]. 
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Restarting the process when it hits one of the x± , it foUows from the strong Markov 
property that we can choose — ^^yg, which concludes the proof by setting 
Ap = Taa,x{A'p}. □ 

3. Convergence of the laws 

As usual in the theory of homogenization, we do not show directly that the 
processes X'^ converge to a limit, but we first introducing a compensator g: M — )■ M 
that 'kills' the strong drift of the rescaled process and consider instead the family 
of processes 

Y^t)^X^t) + eg[^) . (3.1) 

Since we will choose 17 to be a bounded function, the weak convergence in the space 
of continuous functions of the laws of X'^ to some limiting process is equivalent to 
that of the Y". 

In order to construct g, let Ci denote the generator of the diffusion with drift 
bi, that is Ci = \d1 + bi{x)dx, and denote by fii{dx) = exp{2Vi{x)) dx the 
corresponding invariant probability measure on [0,1]. We then denote by gi the 
unique smooth function solving 

C-igi = -hi , / gi{x) ^.i{dx) = . (3.2) 
Jo 

Since h is assumed to be centred on either side of the interface, such a function exists 
(and is unique) by the Fredholm alternative. We now choose any smooth function 
(7: M — > M such that g{x) = g-{x + rj) for x G (—00, —77) and g{x) = g+{x — rf) for 
X € (?7,oo), with a smooth joining region in between. 

The main ingredient in our proof of convergence will be |FW93[ Theorem 4.1], 
which is used in conjunction with the previous tightness result to identify the weak 
limit points of the family of probability distributions as the solutions to a martingale 
problem. The aim of this section is to explain how to fit our problem into the 
framework of }FW93| and to verify the assumptions of their main convergence 
theorem. 

Before we proceed, let us recall what is understood by the "martingale problem" 
corresponding to some operator A (see for example |EK86j ). and let us try to guess 
what the generator A for the limiting process is expected to be. Let X be a Polish 
(i.e. complete separable metric) space; C [ 0, c») , the space of all continuous functions 
on [0, 00) with values in X. For any subset / C [0, 00), denote by ^/ the cr-algebra 
of subsets of C[0, 00) generated by the sets {x € C[0,oo) : x{s) € B}, where s e / 
and i? C X is an arbitrary Borel set. We also denote by C{X) the space of all 
continuous real- valued functions on X. 

Let A be a linear operator on C{X), defined on a subspace ^{A) C C{X). We 
will say that a probability measure P, on (C[ 0, 00), =^[0, 00)) ) is a solution to the 
martingale problem corresponding to A^ starting from a point G X, if 

P{a; : a;(0) = xo} = 1 (3.3) 

and, for any / G ^{A), the random function defined on the probability space 
(C[0,oo),^[o,oo),P) by 

f{x{t)) - f (Af) {x{s)) ds, t G [0, 00) , (3.4) 
Jo 
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is a martingale with respect to the filtration {=^[o,t] }j>q- 

What do we expect the operator A to be given by in our case? On either side 
of the interface, we argued in the introduction that the limiting process should be 
given by Brownian motion, scaled by factors C± respectively. Therefore, one would 
expect A to be given by 

(Af){x) - <^ ? : " (3.5) 
1 \Cldlf{x) otherwise, ^ ' 

and the domain ^^{A) to contain functions that are away from the origin. This 
however does not take into account for the "skewing" , which should be encoded in 
the behaviour of functions in S!{A) at the origin. 

Since the limiting process spends zero time at the origin (the invariant measure 
is continuous with respect to Lebesgue measure), it was shown in | Lej06| that the 
possible behaviours at the origin are given by matching conditions for the first 
derivatives of functions belonging to 2i{A). We know from the introduction that 
the invariant measure of the limiting process is proportional to Lebesgue measure 
on either side of the origin, with proportionality constants A± . We should therefore 
have the identity 

A_ / Af{x)dx + \+ / Af{x)dx = Q, 



for every function / e 2!{A). Using (3.51, we thus obtain 

f'{x)dx = 



/O nc 
f'{x)dx + X+Cl 
-oo Jo 



Integrating by parts, this yields (for say compactly supported test functions /) the 
condition 

= \+ClfiO+) . (3.6) 

This is exactly the general form of a generator produced by Theorem 4.1 in Frei- 
dlin and Wentzell |FW93j . a differential operator on the regions away from some 
distinguished points termed nodes, combined with a restriction on the ratios of the 
limits of the derivatives at this point. 

The main theorem of this section that is also very closely related to the main 
theorem of the article is as follows: 

Theorem 3.1. Let C± be given by 

Cl = f\l + g±(xyf fi±{dx) , (3.7) 

"'0 



where g± and fi± are as in (3.2). Let A be given by \3.^ and let S>{A) be the set 
of continuous functions / : M — >■ M, vanishing at infinity, that are away from 



and that satisfy the condition (3.6) at the origin. 

Then, every limit point of the family of processes Y,f. is solution to the martingale 
problem corresponding to A. 

As already mentioned, our main ingredient is jFWQS) Theorem 4.1] applied to the 



sequence of processes as defined in (3.1 1. For completeness, we give a simplified 



statement of this result adapted to the situation at hand: 
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Theorem 3.2 (Freidlin & Wentzell). Let Ci, i = ±, he elliptic second order differ- 
ential operators with smooth coefficients on li, — [0,oo), /_ = (— oo,0], and let 
be a family of real-valued processes satisfying the strong Markov property. For 
some fixed fj > 0, let be the first hitting time of the set {—ef},ef}) by . 

Assume that there exists a function k: M_|- — > with lime_j.o ^(s) = such that, 
for any function f G C^{Ii) and for any A > 0, one has the bounq}\ 



^'f{Y^T^))^f{y) + 







Y^t))-Cj{Y^t)))dt 



^0{k{e)) 
(3.8) 

as e — >■ 0, uniformly with respect to y £ li. Assume furthermore that there exists a 
function 6: — > with Iuh^^q S{e) — and lini£_j.o (5(£)/A;(e) — >■ oo such that, 
for any A > 0, 

/•oo 



E, 







e-^h^,s,s){Y'{t))dt 







(3.9) 



as £ — ^ 0, uniformly over all y £ R. Finally, writing for the first hitting time of 
the set {—CO, —S) U {S, oo) by Y'^ , assume that there exist Pi > with p_ +p+ = 1 
such that 

¥y[Y'{a')Gl,]^p,, ie {+,-}, (3.10) 

uniformly for y € (—£??, sfj). 

Let now A be the operator defined by Af{x) — Cif{x) for x £ li with domain 
S!{A) consisting of functions f such that /I/; G C^(/i) and such that the 'matching 
condition' p_|_/'(0+) — P-f'{0^) holds. Then for any fixed > 0, any A > 0, and 
any f G S>{A), the bound 



ess sup 



Xf{Y%t))-Af{Y^t)) 



^1 



dt^e-^''>f{Y'{t„)) 





(3.11) 



holds as £ — > 0, uniformly for all y G 



Remark 3.3. The version of Theorem 3.2 stated in |FW93| does actuaUy treat more 
general diffusions on graphs, but assumes that the edges of the graph are finite. This 
is not really a restriction, since J+ is in bijection with [0, 1) (and similarly for /_) 
and we can simply add non-reachable vertices at ±1 to turn our process into a 
process on a finite graph. 



Remark 3.4. As can be seen by combining (1.3) and (3.6|, the probabilities p± 



appearing in the statement of Theorem 1 3. 2 1 are not quite the same in general as the 
probabilities p± = {p, l—p} appearing in the construction of skew Brownian motion 



in the introduction. The relation between them is given by 



The 



{l-p)C- - p- 

reason is that p± give the respective probabilities of hitting two points located at a 
fixed distance from the 'interface', whereas the non-trivial scaling of the Brownian 
bridges on either side of the interface means that p± give the probabilities of hitting 
two points whose distances from the interface have the ratio C-\-/C-. 



Most of the remainder of this section is devoted to the fact that: 



We denote by CJf the space of smooth functions that vanish at infinity, together with all of 
their derivatives. 
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Proposition 3.5 . T he family of processes given by (3.1) satisfies the assump- 
tions of Theorem 3.2 with C± = ^C±d^ and p± defined by the relations 



P- = 1 , 



P- 



and X± as in {1.4 )■ 
This yields the 



Proof of Theorem 3.1 Before we start, let us remark that the initial condition y 
for the corrected process and the initial condition x for the original process 
are related hy y = x + £g(x/e). 

Note also that, thanks to the identity e~^'^F[s) ds = Ae"'*'* rF{s) ds dt 
valid for any bounded measurable function F, the left hand side in ( |3.11 1 can be 
written as 

ft 



A(e) 



to 



Xe-^%(fiY^t))-f{Y%to))- / Af(Y-{s))ds 

^ Jto 

-^%(g/(r^^o,^)|=^[o.to])rf^ ■ 



[0,to] 



dt 



Ae 



We have already established the weak precompactness of the family {P^.e, £ > 0} 
in the space C(M+,E). The uniformity in x of the convergence of A(£) to then 
implies that for any n, any < ti < ■ ■ ■ < t^ < tg, and any bounded measurable 
function G{xi, . . . ,Xn), Xi € M, 



E. 



,(G(y^(ii),...,y^(t„))-^ Ae-^*5/(y^^o,^)d^)| < sup|G| • A(e) . (3.12) 

If we furthermore assume that G is continuous, then the expression inside the 
expectation is a continuous function on C(M_|_,M), so that any accumulation point 
satisfies 

Ae-^*E(G(X°(ti), . . . , X"(^„))g^(X°, io, dt = Q. (3.13) 

Since the integrand is a continuous function of t and a continuous function is deter- 
mined uniquely by its Laplace transform, this implies that E(G(. . .)Qf{X^^ Iq, t)) = 
for all n and < ti < • • • < t„ < so that in particular the random function 
f{X'^{t)) — Af {X^ (s)) ds is indeed a martingale in the filtration generated by 
the process X^. 

Since the laws of the starting points of X^ are all equal to 6^ by construction, 
we conclude that the law of X^ is indeed a solution of the martingale problem 
corresponding to A, starting from xq. □ 



Proof of Proposition \3.5\ The proofs of ( 3.8 ), ( 3.9 ) and ( 3.10 1 will be given as three 
separate propositions. 



Proposition 3.6. There exists fj > such that the process Y'^(t) satisfies 
with k{e) — £, that is, 



E, 



-A 



^'f{Y^{T-))~f{x) + 



,-At 



\f{Y^-{t)) - \clf"{Y^{t)))dt 



= e) 
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for every function f G C^(M+), uniformly in y ^ [rye, oo), and similarly for the left 
side of the interface. 

Proof. The treatment of both sides of the interface is identical, so we restrict our- 
selves to M+. As before, the initial condition y for the corrected process Y"^ and the 
initial condition x for the original process are related by y = a; + eg{x/e). Note 
that one has g'±{x) 7^—1 for any x since otherwise, by uniqueness of the solutions 
to the ODE g" = —2(1 + g')b, this would entail that g±{x) = —1 over the whole 
interval [0, 1], in contradiction with the periodic boundary conditions. 

By possibly making rj slightly larger, we can (and will from now on) therefore 
assume that g'{x) > — 1 uniformly over a; G K, so that the correspondence x ^ y 
is a bijection. Since g is bounded, this shows that one can find 77 > so that 
y ^ [—£77, £77] implies that x ^ [—£?], £77]. In particular, fixing such a value for 77 
from now on, we see that the drift vanishes in the SDE satisfied by , provided 
that we consider the process only up to time t'^. 

Using the integration by parts formula and Ito's formula for each we get 







f{Y.^{T-))=f{x)+ I e-^^{l + g'^{e-'Xl{s)))f'{Y,^{s))dB, (3.14) 

ds 







M{Y^{s)) + \{l + g'^{e-'Xl{s))ff"{Y^{s)) 



Since the expectation of the stochastic integral vanishes (both /' and g'j^ are 
uniformly bounded), all that remains to be shown is that the last term in the above 
equation converges at rate e to the same term with (1 + g')^ replaced by C\. 

This will be a consequence of the following result (variants of which are quite 
standard in the theory of periodic homogenization), which considers the fully pe- 
riodic case. It is sufficient to consider this case in the situation at hand since we 
restrict ourselves to times before r^, so that the process does not 'see' the interface. 

Lemma 3.7. Let b: M" — > M" be smooth and periodic with fundamental domain 
A C M", and denote by fi the (unique) probability measure on A invariant for the 
SDE 

dX{t) = b{X{t))dt + dBt , X{0) = x, (3.15) 

where B is a standard d-dimensional Wiener process. Assume furthermore that 
J^b{x) ^(dx) ~ 0. (This condition will be referred to in the sequel as b being 
centred.) 

Let h ; M.^ M. be any smooth function that is periodic with fundamental domain 
A and centred. Let furthermore X'^(t) — eXit/e^) and let be a family of (possibly 
infinite) stopping times with respect to the natural filtration of B. Then, for every 
F e Co(K'',K) there exists C > independent of such that the bound 



e~^^F{X'{s))h(^^^^ds 



< Ce 



holds for any e £ (0, 1], uniformly in x. 



Proof. Denote by £ = + (fo(x), V) the generator of (3.15) and let g be the 



unique periodic centred solution to Cg = h. (Such a solution exists by the Fredholm 
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alternative.) Applying Ito's formula to the process 



we obtain the identity 



Jl -At 

e e 



F{X^{s))h 



X^{s) 



-Xe-^^F{X{s))g[ 



ds 



X^{s) 



^2 / e~^'AF{X%s))g[ 



VF{X%s))g 
fX^{s) 



ds 



e 



ds 



-\s 



VF{X^s))g[ 



-As 



F{X%s))Wg 



X^is) 



ds 
dBs 
dB, 



-As 



\IF 



(rw).v.(^^) 



ds 



The claim then follows by taking expectations and noting that all the functions of 
X^ appearing in the various terms are uniformly bounded. □ 



Returning to the proof of Proposition 3.6 we first note that f"{Yy) — f"{X^) + 
0{e), so that we can replace f"{Yy) by f (X^) in (3.14), up to errors of 0{e). 
Applying Lemma 3.7 with F = f" and h = {1 + g'_^)^ — C^, the claim then follows 
at once. □ 



Proposition 3.8. The convergence 



E 



-Ati 







(3.16) 



takes place as e — > 0, uniformly in the initial point y e M. In particular, ( 3.9) holds 
with 6{e) = \j£. 

Proof. The main idea is to first perform a time-change that turns the diffusion 
coefficient of into 1 and to then compare the resulting process to the process 
which is the solution to 

dV^ = bl^iV) dt + dBt , (3.17) 
where the drift by is given by 

{-^ forO<x<77e, 
^ for -?}£ < X < 0, (3.18) 

otherwise, 

for Cv and rj some positive constants independent of e to be determined below. An 
explicit resolvent equation then allows one to show that (3.16) with Y replaced by 
V tends to zero as e — ?> 0, uniformly in the initial point. 
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First, let us start with the time change. As in the proof of Proposition |3.6[ 
we assume that g is chosen in such a way that g' is bounded away (from below) 
from —1, so that there exists a constant > such that g'{x) > ^ — 1 for every 
cc € M. In order to turn the diffusion coefficient of Y"^ into 1, we use the time change 
associated with the quadratic variation of the Y^ , 



{Y^Y^m 



ds 



thus setting 



Defining the function b 
equation 



CI ^ ini{t' > : {Y',Y'){t') > t} . 

b + Cg, the process Z^{t) = y^(Cf ) then satisfies the 



dB, 



Note that the time-change was defined precisely in such a way that the second term 
in this expression is equal to some Brownian motion B'^ (t) . Inserting the expression 
for the time change, the first term can be rewritten as 



rh(^)ds=rh( 

J„ e \ e J 7o £ ^ 



1 



9 



ds 



It follows that the drift term is non-zero only when the time-changed process occu- 
pies the region (— £77 — e||(7||oo, £?? + e||5l|oo) just as for the non time-changed process. 
To summarise, there exists a drift b bounded uniformly by ^ for some constant 
Cv > and vanishing outside of [—fie, fje) for fj = rj+ Uglloo, as well as a Brownian 
motion B^, so that the process satisfies the SDE 



rfZ^ = fe(Z^) dt + dBI 



zm = y 



(3.19) 



Now, look at how the time change affects the expression (3.16), where we set 

{Y^t))dt 



T(t))dt 



e ice 



< sup(l 



lG^{Z%t))[l + g 



dt 



e ice 



'Z^{t))dt 



Hence if it can be shown that 



e-^^lG^-{Z^{t)) dt 







(3.20) 



uniformly in the initial point x for the underlying process X|, as e — >■ 0, then our 
claim follows. The idea is to bound (3.20) by the 'worst-case scenario' obtained by 
replacing the process Z'^ by the process described in (3.17). 

One technical problem that arises is that it is tricky to get pathwise control on 
the behaviour of V due to the discontinuity of its drift. We therefore first compare 
Z'^ with the process /7| solution to 



dW = ¥jj{U') dt + dBlj{t) , U'^iO) = y , 



(3.21) 
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where Bfj is a Brownian motion to be determined and bfj is the Lipschitz continuous 
odd function defined on the positive real numbers by 



Cv 



-4((3 + ^)e- 



x) 







for |x| < e, 

for e < X < (2 + 77)2, 

for {2 + f])e < X < {3 

otherwise. 



■V)e, 



The SDE (3.21) satisfies pathwise uniqueness, which is why we are using it as 
an intermediary between Z'^ and V^. Now, what we are going to do is, given a 



reahsation of the Brownian motion driving in (3.191, to choose the Brownian 



motion Bfj driving by changing the sign of the increments in such a way that 
the absolute value of is always less than or equal to \Z^\ + 2e. By pathwise 



uniqueness, we are indeed free to choose the Brownian motion in (3.21). The 
choice of the Brownian motion is the content of the following lemma: 

Lemma 3.9. For every initial condition x, there exists a map i-> Bfj that pre- 
serves Wiener measure and such that \U^\ < \Z^\+2s almost surely. In particular, 
it follows that ( 3.2(Jj^ is bounded by 



Ex 



(3.22) 



where S' = 6 + 2e. 



Proof. The construction works in the following way. Consider first the processes 
driven by the same realisation B^ and define a stopping time tq by tq = mf{t > 
: |C/|(t)| ~ \Z^{t)\ + 2e}. This stopping time is strictly positive and one has 
|C^:^('''o)| > 2£. For times after tq, we determine Bfj by 

Bfj{t) = Bfjiro) + szgn{U^iTo)) f szgn{Z'^{s)) dB'{s) , 



and we introduce the stopping time ai = iid{t > tq : \U^\ — e}. Since by construc- 
tion C/^ does not change sign between tq and <Ji , it then follows from the Ito-Tanaka 
formula that up to cti one has 

d\U'\ = 6^(|C/"|) dt + sign{Z') dB^{t) , 

d\Z^ = sign{Z^y¥ {Z') dt + sign{Z^) dB%t) + dL{t) , 

for some local time term L. Since the local time term always yields positive con- 
tributions and since it follows from the definition that bfj{u) < sign(z)b'^ (z) for 
|u| > e and \u\ < \z\ + 2e, we can apply a simple comparison result for SDEs to 
conclude that the inequality \U^\ < \Z^ \ +2e holds almost surely between times tq 
and do- 

We then drive again both processes by the same noise and define as before ti by 
n = mi{t > ai : |f/^(t)| — \Z^{t)\ -|-2e}. Note that ti > tri almost surely since one 
has |t/^(CTo)| < l^x(o'o)| + £■ We then apply the previous construction iteratively, 
so that, setting an = 0, we have constructed Bfj by 



Bf 



n=0 71=0 
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Since the process U' has finite quadratic variation and has to move by at least e 
between any two successive stopping times, our sequence of stopping times does 
converge to infinity, so that Bfj[t) is indeed a Brownian motion with the required 
property. □ 

In our next step, we compare the process that we just constructed with the 
process defined in (3.17), where we set 7) = 5 + i). Since the drift coefficient 



is bounded, it follows from an application of Girsanov's theorem like in [ RY9H 
Corollary IX. 1.12] that this SDE has a solution for some Brownian motion By, say. 



We now fix By and use it to construct a Brownian motion Bfj driving (3.21 ) in 
such a way that the absolute value of always stays less than | 



2e: 



Lemma 3.10. 

such that \V^\ 
hounded by 



There exists a map By 1— >■ Bfj that preserves Wiener measure and 
— I^^kI + 2£ for all times almost surely. In particular, (3.22) is 



e-^*l 



(-5", 5") 



{V^{t))dt 



(3.23) 



with S" = S' + 2e. 



Proof. The argument is virtually identical to that of Lemma |3.9[ so we do not 
reproduce it here. □ 

It now remains to show that: 



Lemma 3.11. The expression (3.23) converges to uniformly in the initial point 
as e — > 0. 

Proof. We write e for 5e + er] and 6 for 5 + Ae for ease of notation, but this has no 
bearing on the rates of convergence of the aforementioned quantities and hence on 
the calculation. We have the identity 



E, 



e-^'li-s,S){V^{t))dt 



-At 



PHI 



) [x] dt 



^{\-C\.)-\^sMx) (3.24) 

by the resolvent equation (see for example |EK86|, Chapter 1]), where Cy is the 
generator of the Markov semigroup P^ associated to . 

We now proceed to computing this expression explicitly in order to show that 
its supremum tends to zero uniformly in x. In order to keep notations simple, we 
assume for the remainder of this proof that 77 — Cy — 1, which can always be 
achieved by rescaling space and redefining e. In this case, the solution f{x) = 
(A — £y)^^l(_5^5)(a;) to the resolvent equation can be assembled piecewise on the 
intervals {—oo,S), (— (5, — e), (—£,0), (0,£), {e,S) and {S,oo) by making sure that 
it is at each junction. Owing to the symmetry of the problem, the function / 
will be an even function of x, hence we only have to analyze it on one side of the 
origin. 

The general solution on each interval can be written as 



/(^) = 



iJoe""^ for x> S, 

i +ylie^^ + Sie-^^ fore<a;<5, 
i + e'^A2e'^'^'' + B2e-^''= for x < e, 
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where 



71 



72 



2A 



2A 



0{s) 



The reason for the somewhat strange choice of adding an expUcit factor in front of 
A2 is justified a posteriori by noting that with this scahng, the matching conditions 
at £ and 5 (as weU as the fact that the derivative should vanish at the origin) yield 
the following linear system: 



M 



fBo\ 




( ^ 












2 


-A\ 


Ai 












-1 


-1 





1 


Bi 







M = 





-1 


1 








A2 




1 

A 




-1 


1 


1 


















1 


-1 





0/ 



0(5) 



To lowest order in e and 5, this can easily be solved exactly, yielding 
(50,^1,51,^2,52) = -^(0,1,1,A,2) +0(<5) . 

Inserting this into the expression for / shows that sup^g^ |/(x)| — 0{5) — 0{^/e), 
thus completing the proof. □ 



sit ion 



With the lemma regarding the resolvent calculation above, the proof of Propo- 

□ 



3.8 is complete. 



We finally show that 

Proposition 3.12. For every c > 0, the exit probabilities from the interval (—5, 5) 
satisfy the bound 

uniformly for x G [— c£, ce] . 

Proof. For the proof of this result, it turns out to be simpler to consider the original 
process X^{t). 

Whenever Y^{t) exits the set {—S,S), due to the deterministic relationship be- 
tween the processes, X^{t) exits a set {—S',S"), where S' and 5" are contained in 
the interval {S — e||(7||oo,<5 + e||g||oo)- Therefore, we just look at the exit of X^{t) 
from an interval of this form as the computations are much easier to carry out. 
This is due to the simpler form of the scale function for X'^{t) compared with that 
of Y^it). It follows from |RY91[ Exercise VII.3.20] that the scale function of the 
diffusion on R with generator £ = ^ cr^ (x) 



with generator C 

S{X) : 



^(^) S is gi'^sn by. 



exp 



2b{z)(j^\z)dz) dy 



where c is an arbitrary point in M. Recall that the scale function of a real-valued 
process is a continuous, strictly increasing function such that for any a < x < 6 in 
the set where the Markov process takes its values, one has 

s{x) — s{a) 



s[b) — s{a) ' 
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where Ta , are the first hittmg times of the points a and b respectiveiy. For (t) 
we have that a = 1 and the drift is equal to ^ b{x/e). We are from now on going 

to use the notation qe{y) = (l +eg{-/£)) ^(y) for the transformation that allows 
recovery of X"^ from . We also denote by the first hitting time of the point 
a G M by the process X^ . We also use the shorthand notation 



F^{u) = exp(^-2 j - b{z/e) dz 



so that the scale function for X^ is given by s{z) = Fe(y) dy. 

With this notation at hand, we have, for x G {—S,S), that, denoting the escape 
time of Y'^{t) from {—S, S) by cr^ , 

/J'-^-) dy - j;'' FM dy Ilf^ F^jy) dy 

j;'"FMdy-j;''FMdy f:'s.FMdy ' 

Noting that F^ has the scaling property F^{u) — Fi{u/e), we thus obtain the 
identity 

. rf^Fi{y/e)dy j\,F,{y)dy + Oil) 

Cs'F^(^yl^)dy J%F,{y)dy + 0{l) 

where we used the fact that qe{x) — 0{e) and Fi is uniformly bounded, due to the 
fact that the functions b± are centred by assumption. 

Note now that the effective diffusion coefficients C± can alternatively be ex- 
pressed as [PS081 Sec 13.6] 



i?+l rV+i 

exp(-2V^(u)) du / exp(2y(M)) du 



and similarly for C_. Therefore, since Fi is periodic away from [— ?7,7?], it follows 
immediately from the definitions of \± and C± that 



F^{y)dy^-^N + 0{1), / F^{y) dy = N + O {I) . 

Since 6 = ^/e, combining these bounds with ( |3.25 ) imphes that fxiY^i^cr^) e /+) = 
C\X+/{C'i\+ + C^A_) -I- C(v/e), from which the requested bound follows. □ 

Combining Propositions |3.6[|3.8[ and |3.12| completes the proof of Proposition |33] 

□ 



4. Uniqueness and characterisation of the martingale problem 

To conclude this article, we show that solution to the martingale problem cor- 
responding to A is unique and is indeed given by the variant of skew Brownian 
motion constructed in the introduction. 

The skew Brownian motion Bp of 'skewness' parameter p is known to have gen- 

erator Cp = ^ on the set of functions that are continuous, twice continuously 
differentiable except at the origin where we have pf'{0~^) = (1 — p)/'(0~). (See for 
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example the review article |Lej06| .) The process Bc^^p constructed in the intro- 
duction is given by 

Bc,At)-Gmt))^ = othe^wL 

Since G is a continuous bijection, this is again a strong Markov process and it has 
generator given by Af = (£p(/ o G)) o G"! with 2i{A) ^ {f : / o G € ^(£p)}. 



Using the relation (1.3), it is now a straightforward calculation to show that ^^{A) 



consists precisely of those functions satisfying the derivative condition in (3.6), 
hence we will have 

P+ ^ C+P 
p_ G_(l-p)- 

This provides a clue as to how to show uniqueness easily. Let Y denote any solution 
of the martingale problem corresponding to A acting on &{A), Y also represents 
any possible limit point of the family as e — ^ 0. Defining g like G, but with 
constants p/C+ and (1 — p)/C- instead of G_|_ and G_, we note that g satisfies the 
derivative condition at imposed for elements of &{A). Since g doesn't vanish at 
infinity, we have g (fi ^{A), but we can approximate g by a sequence gn € 3>{A) 
such that g = gn on [— n,n]. Indeed, given i > 0, the probability of escaping from 
[— n, n] before time t tends to zero as n — >■ oo uniformly for e G (0, 1] by tightness. 
As a consequence, the process V = g{Y) satisfies the SDE 



Jo 



where A± are some constants. It is known |BC051 Theorem 2.1] that this equation 
has a pathwise unique solution which in particular implies uniqueness in law. Since 
g is an invertible map, this immediately implies that Y is unique in law and one 
can check that Y is indeed the variant of skew Brownian motion constructed in the 
introduction. 
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